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Abstract. It is thought that a stochastic background of gravitational waves was 
produced during the formation of the universe. A great deal could be learned by 
measuring this Cosmic Gravitational- wave Background (CGB), but detecting the 
CGB presents a significant technological challenge. The signal strength is expected 
to be extremely weak, and there will be competition from unresolved astrophysical 
foregrounds such as white dwarf binaries. Our goal is to identify the most promising 
approach to detect the CGB. We study the sensitivities that can be reached using both 
individual, and cross-correlated pairs of space based interferometers. Our main result 
is a general, coordinate free formalism for calculating the detector response that applies 
to arbitrary detector configurations. We use this general formalism to identify some 
promising designs for a GrAvitational Background Interferometer (GABI) mission. 
Our conclusion is that detecting the CGB is not out of reach. 
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1. Introduction 



Before embarking on a quest to discover the Cosmic Gravitational-wave Background 
(CGB), it is worth reflecting on how much has been learned from its electromagnetic 
analog, the Cosmic Microwave Background (CMB). The recent Boomerang|I| and 
Maxima[0| experiments have furnished detailed pictures of the universe some 300,000 
years after the Big Bang. These measurements of the CMB anisotropics have been used 
to infer that the visible universe is, to a good approximation, spatially flat. The earlier 
COBE-DMR0 measurements of the CMB anisotropy flxed the scale for the density 
perturbations that formed the large scale structure we see today, and showed that the 
density perturbations have an almost scale free spectrum. With all this excitement 
over the anisotropy measurements it is easy to forget how much the CMB taught us 
before the anisotropics were detected. The initial observation that the CMB is highly 
isotropic is one of the major observational triumphs of the Big Bang theory. Equally 
important was the discovery by COBE-FIRAS[Q that the CMB has an exquisite black 
body spectrum with a temperature of 2.728 Kelvin. 

The lesson that we take from the CMB is that while it would be wonderful to have 
a COBE style map of the early universe in gravitational waves, a great deal can be 
learned by detecting the isotropic component. The main focus of this work is on the 
flxing the amplitude of the CGB at some frequency. Some attention will also be given 
to the prospects of measuring the the energy spectrum. Detecting anisotropics in the 
CGB is a considerably harder problem that we address elsewhere 0. 
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Figure 1. Possible sources of a stochastic gravitational wave background plotted 
against the sensitivity curve for the LISA mission. 



There are two major obstacles that stand in the way of detecting the CGB. The 
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first is the extreme weakness of tlie signal, and tlie second is tlie competing stocfiastic 
background produced by astropfiysical sources. Figure 1 sliows various predictions for 
the CGB energy spectra in several early universe scenarios including standard inflation, 
a particular M-theory model^ and quintessence based inflation|]^. Also shown is the 
sensitivity curve for the proposed Laser Interferometer Space Antenna (LISA) , and a 
compilation of possible extragalactic astrophysical foregrounds taken from the work of 
Schneider et a/.[^. The spectrum is expressed in terms of f2gw(/), the energy density in 
gravitational waves (in units of the critical density) per logarithmic frequency interval, 
multiplied by /iq, where is the Hubble constant in units of 100 km s~^ Mpc~^. The 
predictions for the CGB power spectra are fairly optimistic, and some of them come 
close to saturating existing indirect bounds on Qg^{f) - see the review by Maggiore for 
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With the exception of the M-theory motivated model, we see that the major 
obstacle to detecting the CGB is not the sensitivity of our detectors, but competition 
from astrophysical foregrounds. The astrophysical foregrounds are produced by the 
combination of a great number of weak sources that add together to form a significant 
stochastic signal. Based on current estimates, it is thought that the contribution from 
astrophysical sources falls off fairly rapidly below 1 yuHz. For that reason, we will 
focus our attention on detecting the CGB in the sub-micro Hertz regime. To put this 
number into perspective, an interferometer built using LISA technology would need 
spacecraft separated by one-tenth of a light year to reach a peak sensitivity at 1 /xHz! 
A more practical approach is to use two smaller interferometers that lack the intrinsic 
sensitivity to detect the CGB on their own, but together are able to reach the necessary 
sensitivity. The idea is to cross correlate the outputs from the two interferometers and 
integrate over some long observation time T. Since the noise in the two interferometers is 
uncorrelated while the signal is correlated, the signal to noise ratio will steadily improve 
as the observation time is increased. 

In principle, it is possible to detect an arbitrarily weak signal by observing for an 
arbitrarily long time. In practice, the prospects are not so rosy as the sensitivity of the 
detector pair only improves on the sensitivity of a single detector as (TA/)^/^, where 
A/ is the frequency bandpass. The frequency bandpass is typically taken to equal 
the central observing frequency. Suppose we try and look for the CGB at a frequency 
of / ~ A/ = 10^ Hz using the two LIGO detectors. Over an observation time of 
T = 1 year, the sensitivity is improved by a factor of ~ 2 x 10^ by cross correlating 
the Washington and Louisiana detectors. For a pair of LISA detectors with a bandpass 
of A/ ^ 10"^ Hz, cross correlating for a year would result in a 500 fold improvement 
in sensitivity. However, for detectors operating in the /iHz range, cross correlating two 
detectors for one year only improves on the sensitivity of a single detector by a factor 
of ~ 6. So why use two interferometers when one will do? The reason is simple: with 
a single interferometer it is not possible to tell the difference between the CGB and 
instrument noise since both are, to a good approximation, stationary Gaussian random 
processes. Moreover, it is impossible to shield a detector from gravitational waves in 
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order to establish the noise floor. With two or more independent interferometers the 
cross correlated signal can be used to seperate the stochastic bacground from instrument 
noise. The LISA system has three arms, and the signal from these arms can be used 
to form three different Michelson interferometers. However, these interferometers share 
spacecraft, and thus common sources of noise, so nothing is gained by cross correlating 
the outputs. However, Tinto et a/.[rT| have recently suggested another way of combining 
the signals, known as a Sagnac system, that creates an interferometer that is fairly 
insensitive to gravitational waves. Since the response of this "bad" interferometer is 
noise dominated, its output can be used to estimate the noise floor for the standard 
interferometer configuration. Thus, a single LISA type observatory might be able to 
discriminate between the CGB and instrument noise. 

We will show that cross correlating two LISA interferometers results in a signal to 
noise ratio of 

SNR^1.4m"7^) (1) 

Vyear / V I 

for a scale invariant stochastic gravitational wave background. This represents a five 
hundred fold improvement over the sensitivity of a single LISA interferometer. If it were 
not for the astrophysical foregrounds that are thought to dominate the CGB signal in 
the LISA band, a pair of LISA detectors would stand a good chance of detecting the 
CGB. In an effort to avoid being swamped by astrophysical sources, we focus our efforts 
on detecting the CGB below 1 yuHz. We show that a pair of interferometers can reach 
a signal to noise of 

for a scale invariant CGB spectrum. The above result is scaled against a possible LISA 
follow-on mission (LISA II) that calls for two identical interferometers, comprising six 
equally spaced spacecraft that form two equilateral triangles overlayed in a star pattern. 
The constellations follow circular orbits about the Sun in a common plane at a radius 
of 1 AU, so that each interferometer arm is L = AU in length. The acceleration 
noise, 5a, is scaled against a value that improves on the LISA specifications by one 
order of magnitude. The LISA II design could detect a stochastic background at 99% 
confidence with one year of observations if f2gw^o ~ 10~^^. Ideally we would like to 
reach a sensitivity of at least Vt^^K^ ~ lO^^*' in order to detect inflationary spectra with 
a mild negative tilt. It is difficult to achieve this with spacecraft orbiting at 1 AU as 
thermal fluctuations due to solar heating of the spacecraft make it hard to reduce the 
acceleration noise below 5a ~ 10~^^ms~^. The best strategy is to increase the size of 
the orbit, i?, as this increases the size of the interferometer, L = \/3R, and decreases 
the thermal noise by R~^. 
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2. Outline 



We begin by deriving the response of a pair of cross correlated space-based laser 
interferometers to a stochastic background of gravitational waves. Similar calculations 
have been done for ground based detectors [|1^, |13|, [1^ |15|, [I6[, but these omit the high 



frequency transfer functions that play an important role in space based systems. Our 
calculations closely follow those of Allen and Romano |[16[|, and reduce to theirs in the 
low frequency limiti]]. Many of the calculational details that we omit for brevity can 
be found in their paper. Having established the general formalism for cross correlating 
space based interferometers we use our results to identify some promising configurations 
for a GrAvitational Background Interferometer (GABI) mission to detect the CGB. 



3. Detector Response 

We attack the problem of cross correlating two space based interferometers in stages. 
We begin by reviewing how the Doppler tracking of a pair of spacecraft can be used 
to detect gravitational waves. Using this result we derive the response of a two arm 
Michelson interferometer and express the result in a convenient coordinate-free form. 
The response of a single interferometer is then used to find the sensitivity that can be 
achieved by cross correlating two space based interferometers. 



3.1. Single arm Doppler tracking 



Following the treatment of Hellings|jT^, we consider two free spacecraft, one at the origin 



of a coordinate system and the other a distance L away at an angle 9 from the z axis. 
Both spacecraft are at rest in this coordinate system. The spacecraft at the origin sends 
out a series of photons, while a weak plane gravitational wave is passing through space 
in the +z direction. To leading order, the spacetime metric in the transverse-traceless 
gauge is: 

ds^ = —c^dt^ + dz^ + [1 + h cos 2ip)dx'^ + {1 — h cos 2ip)dy'^ — h sin 2ipdxdy , (3) 

where h{t — z/c) is the wave amplitude and ip is the angle between the principal 
polarization vector and the x axis. If we choose our coordinates such that the second 
spacecraft is in the x — z plane, the path of the photons can be parameterized by 

a; = p sin 6^ z = p cos 9 . (4) 

The photon path in the perturbed spacetime is given by ds"^ = 0, or 

c dt ±{1 + cos 2ijsm^ 9) dp. (5) 

II Here high and low frequencies are defined relative to the transfer frequency of the detector = 
c/(2ttL), where L is the length of one interferometer arm. 
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The round trip journey from spacecraft 1 to spacecraft 2 and back again is given by 

£{t2-to) = [^\dt = 2L + -cos2ipsm^e I hlto + -{I - cos e)]dp 
Jto 2 \Jo c 

+ r h[ti + + cos9)]dp] . (6) 

Jo c c J 

Here to is the time of emission, ti is the time of reception at spacecraft 2 and ^2 is the 
time of reception back at the first spacecraft. Using h{q) = hQexp{iujq), the varying 
portion of the round-trip distance is 

6i{t2) = Lh{t2)cos2^sm^e^ |^sinc |^^(l-cos^)j exp (^-i^{3 + cosO)^ 

+sinc + cos0)j exp (-^^(1 + cos^)^^ , (7) 

where / = u/{2n) is the frequency of the gravitational wave and /* = c/{27iL) is the 
transfer frequency. To zeroth order in h, the optical path length is i = 2L. The general, 
coordinate independent version of this expression is given by 

Si 1 ^ ^ 

- = -(u®u):h(^])T(u.^],/), (8) 

where u is a unit vector pointing from the first to the second spacecraft and f2 is a unit 
vector in the direction the gravitational wave is propagating. The colon denotes the 
double contraction a : b = aijlfK The transfer function T is given by 



sine ~ " ■ ^) j \ + VL-Vt) 



+ sine ( ^(1 + u ■ Vt)] exp {-i^{l + u-Q) 



(9) 



2fif j y 2/^ 

This expression for T agrees with the one derived by Schilling|]18[. The gravitational 
wave is described by the tensor 

h(t, x) = h-^{ujt - uoVL ■ x)e+(l], i)) + /i^ (cut - uoVL ■ x)e^ {VL, ^) , (10) 

with polarization tensors 

e+(f2,^) = e+{VL) cos2^-e^(f]) sin 2^ 

e^'iVL,^) = e+{n)sm2^p + e''{n)cos2^p. (11) 
The basis tensors can be written as 

e'''(n) = m®m — h®n 

e^{h) =m®h + h®m (12) 

where m, n and VL are an orthonormal set of unit vectors. The expression in (|^) can be 
recovered from (H) by setting 

h'^ = h, = 0, Q = z, u = X sin ^ + i cos 6*, 



x®x — y^y, and = x®y + y®x. (13) 
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The transfer function T approaches unity for f <^ f^, and falls of as 1// for 
/ ^ /* by virtue of the sine function. For a ground based detector such as LIGO 
the transfer function can be ignored since instrument noise keeps the operational range 
below LIGO's transfer frequency of ~ 10^ Hz (not to mention that LIGO is a Fabry- 
Perot interferometer so our results do not really apply[0]). For a space based detector 
such as LISA, the instrument noise does not rise appreciably at high frequencies, and 
it is the transfer function that limits the high frequency response. The sine function is 
responsible for the wiggly rise in the LISA sensitivity curve at high frequency that can 
be seen in Figure 4. 

3.2. Interferometer response 

We can form an interferometer by introducing a third spacecraft at a distance L from 
the corner spacecraft and differencing the outputs of the two arms0|. The measured 
strain is given by 

= D(f!,/):h(n,/,x,i), (14) 

where 

D(f2, f) = ]- ((u ® u) r(u . fi, /) - (v ® v) r(v ■ fi, /)) (15) 



2 

is the detector response tensor and u and v are unit vectors in the direction of each 
interferometer arm, directed out from the vertex of the interferometer. The above 
expression gives the response of the interferometer to a plane wave of frequency / 
propagating in the VL direction. A general gravitational wave background can be 
expanded in the terms of plane waves: 

/oo r ^ ~ ^ 

df I dVL hij{VL, /,x,t) 
-oo J 

= T. r df fdn fiAif, n)e-2-/*e2'^^^^-''/^(fi) . (16) 

A J— OO J 



A 



Here / dil = Jq^ dcj) sin 9 d9 denotes an all sky integral. The Fourier amplitudes obey 
hA{—f) = h\{f) since the waves have real amplitudes. In the final expression we 
have chosen e"*" and as basis tensors for the decomposition of the two independent 
polarizations. The response of the detector (which we are free to locate at x = 0) to a 
superposition of plane waves is then 

/oo r ^ ^ 

df / dn s(n,/,o,t) 
-co J 

df / dn /i^(/,fi)e-2-^^*D(fi,/) : e^ifl). (17) 

^ More complicated differencing schemes liave to be applied if the arms of the interferometer have 
unequal lengths in order to cancel laser phase noise |19[. 
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From this we infer that the Fourier transform of s{t) is given by 

s{f) = J2[d^ hA{L ^^)D(fi, /) : e^{n) . (18) 

In the low frequency hmit the transfer function T approaches unity and /) only 
depends on the geometry of the detector: 

D = ^(u(g)u-v(g)v). (19) 

To a first approximation, a stochastic gravitational wave background can be taken 
to be isotropic, stationary and unpolarized. It is fully specified by the ensemble averages: 

{hAif,n)) = , {hAif,n)hAif',n')) = o , (20) 

and 

^ 1 pfO O'] 

{h*M,^)hA'{f\n')) = -6{f - fY-^^6AA' S,{f) . (21) 

Here Sh{f) is the spectral density of the stochastic background and the normalization 
is chosen such that 



{h*A{f)hAif)) ^ J dQdQ'{h\{f,Q)hA'{f,n')) 

= \5U -n^AA'SnU)- (22) 

The spectral density has dimension Hz"^ and satisfies Sh{f) = Sh{—f)- It is related to 
^gw(/) by 

47r2 

^^gw(/) = ^/'^.(/). (23) 
The time averaged response]^ of the interferometer, 

(s(t))r = lim - s(t)dt, (24) 
when applied to a stochastic background, is equivalent to the ensemble average 

= {S{t)) 

df / dn {hA{f,n))e-'-''^'-D{Q,f) : e^Q) 

= 0. (25) 

Since the expectation value of sit) vanishes we need to consider higher moments such 
as s^(t). Using equations ([T^) and (|l^) we find 

{s\t)) = / dfSuUmf) , (26) 
where the transfer function 7^(/) is given by 

nf) = / ^ E f)F\^, fT , (27) 

The limit r ^ oo is approximated in practice by observing for a period much fonger than the period 
of the gravitational wave. 
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and F'^{(l^ f) = Y){Q, f) : e^{Q) is the detector response function. In the low frequency 
hmit, / <^ it is easy to show that TZ{f) = 2/5 sin^ /3, where cos /5 = u ■ v is the angle 
between the interferometer arms. 

The response of the interferometer can be expressed in terms of the strain spectral 
density, hs{f), which has units of Hz~^/^ and is defined by 

{s\t)) = / dfhlif) = / dfSnUmf) (28) 



We see that the strain spectral density in the interferometer is related to the spectral 
density of the source by 



hsU) = phUmf) ■ (29) 
The total output of the interferometer is given by the sum of the signal and the noise: 

S{t) = s{t)+n{t). (30) 
Assuming the noise is Gaussian, it can be fully characterized by the expectation values 

(n(/))=0, and (n*(/)n(/')) = i5(/ - /')5„(/) , (31) 

where Sn{f) is the noise spectral density. The total noise power in the interferometer is 
thus 

{n\t))= / dfSn{f)= / dfhlif) (32) 
Jo JO 

where hn{f) is the the strain spectral density due to the noise. Comparing equations 
(p8|) and (|32D, we define the signal to noise ratio at frequency / by 

SNR(/)^|i^l^^fW). (33) 

Sensitivity curves for space-based interferometers typically display the effective strain 
noise 



Snif) ^34^ 



At high frequencies it is the decay of the transfer function TZ{f) leads to a rise in the 
effective noise floor. The actual noise power, S'„(/), does not rise significantly at high 
frequencies for space based systems. 



3.3. The LISA interferometer 

As a concrete example of the general formalism described above we derive the sensitivity 
curve for the LISA interferometer using the design specifications quoted in the LISA 
Pre-Phase A Report [§. Using the coordinate system shown in Figure 2, the unit vectors 
along each arm are given by 




and the gravitational wave is described by 

Q = cos sin 6 ct + sin sin 6y + cos 6 z (36) 
rh = sin (j)x — cos (p y (37) 
fi = cos cos 6 X + sin cos 6y — sin 6' 5 . (38) 

The angle between the interferometer arms is arccos(u ■ \) = j3 = n /3. The various 
ingredients we need to calculate 7^(/) are 

u ■ i7 = sin(0 + 7r/6) sin (39) 
v-n = sin(0 - 7r/6) sin 6 (40) 

and 

1 1 /s 

(u (g) u) : e+ = - sin^ 6* + - cos^ 9 cos^ ^ sin 20 (1 + cos^ 9) . (41) 

1 1 /s 

(v v) : e+ = - sin^ 6* + - cos^ 9 cos^ + ^ sin 20 (1 + cos^ 6') . (42) 

(u® u) : e"" = -cos^sin(20 + 7r/3) . (43) 

(v®v) : e"" = -cos^sin(20-7r/3). (44) 

Before proceeding to find 7^(/) it is instructive to evaluate the detector response 
functions, F^{Q,f), at zero frequency: 

F+(i7) = ^(u (g) u - V (g) v) : e+ 
x/3 

= - ^sin20(l + cos^^), (45) 
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and 

F'^iyi) = ^(u (g) u - V (g) v) : e"" 

= ^ cos 20 cos ^, (46) 

The magnitudes of these response functions are shown in Figure 3. They can be thought 
of as the polarization dependent antenna patterns for the LISA detector appropriate to 
a stochastic background of gravitational wavesQ. 




Figure 3. The magnitudes of the detector response functions F^{n, /) and (il, /) 
in the low frequency Umit. 



Inserting equations (^9D through (^4D into equation (pTD gives the transfer function 
TZ{f) in terms of an integral over the angles 6 and 0. Our explicit expression for LISA's 
transfer function agrees with the one found by Larson, Hiscock and Hellings[^] using 



an alternative approach. We were unable to perform the angular integral to arrive at a 
general form for TZif), but for low frequencies the integrand can be expanded in a series 
to give 

7^^n - 2 sln^ I^l _ 37cos^(/3/2)-10cos^(/j/2)-l / / ^ 
/A/j-^sm/^ i 84cos2(/5/2) [fj 




20 cos'^(/3/2) + 163 cos^(/j/2) - 80 cos^(/j/2) - 6 ( f]' _ , 

2268cos2(/?/2) ' ^ ' ••• ' • ^ ' 

At frequencies above / ~ |/*, the sine function takes over and the transfer function 
falls of as 1//^. Setting /3 = tt/3, we find that a good approximation for the transfer 

* These shapes are not coordinate independent. They depend on how we choose u and v in relation 
to e+ and 
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10 504 \f,) 9072 \f,) 47900160 \f,^ 
16783143 /3/*V 
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[ 126156800 V2/y ' ' ~ ^' 

The coefficient in front of the high frequency term is chosen so that the transfer function 
is continuous at / = 



-16- 




-A -3 -2 -1 

log (f/Hz) 

Figure 4. The effective noise floor for the LISA mission. The solid line was obtained 
numerically while the dashed line is our analytic approximation. 

Now that we have calculated the transfer function, our next task is to estimate the 
the detector noise power, S'„(/). There are many noise contributions discussed in the 
LISA Pre-Phase A Report [^]. The dominant ones are thought to be acceleration noise 
from the inertial sensors and position noise due to laser shot noise. A position noise 
of 5x = 2 X 10^^^ m Hz~^/^ is quoted for each LISA spacecraft. There are two such 
contributions per arm, giving a total of 4 contributions for the interferometer. Since 
the contributions are uncorrelated they add in quadrature to give a total position noise 
of 25x. Dividing this by the optical path length of 2L and squaring gives the position 
noise power: 




An acceleration noise of 5a = 3 x 10 ms ^ Hz is quoted for each inertial 
sensor. This noise acts coherently on the incoming and outgoing signal, for a combined 
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acceleration noise of 25a per spacecraft. There are four such contributions in the 
interferometer that add in quadrature for a total acceleration noise of A5a. Dividing 
this by the square of the angular frequency of the gravitational wave yields the effective 
position noise due to spurious accelerations. Dividing the effective position noise by the 
optical path length and squaring gives the acceleration noise power: 

Adding together the acceleration and position noise power gives the total noise power. 
Using the nominal values for the LISA mission [§] we obtain 

S„(/) = 1^9.24 X 10-^° (^)' + X 10""' j Hz-i . (51) 

Inserting the above expression for Sn{f), along with the approximate expression for the 
transfer function (|4^), into equation (0) yields a useful analytic approximation for the 
effective strain noise, hes{f), in the LISA interferometer. The analytic approximation 
is compared to the full numerical result in Figure 4. 

LISA reaches a peak sensitivity in the frequency range 3 x 10^'^ —>■ 10~^ Hz. Using 
equations ( p3D and (|3^) we have 

S,{f) = hlM)Smf)- (52) 
For a SNR of 2, this translates into a sensitivity of 

hMf) = 6.2 X 10- (%^)"' (^)'" Hz-V^ (53) 

Thus, a stochastic background with Qg^{f)hl > 7 x 10^^^ should dominate LISA's 
instrument noise for frequencies near 3 mHz. The difficulty would come in deciding if the 
interferometer response was due to instrument noise or a stochastic background, as both 
are Gaussian random processes. One way to be sure is to fly two LISA interferometers 
and cross correlate their outputs. 

4. Cross correlating two interferometers 

Having demonstrated that our formalism recovers the standard results for a single 
interferometer, we now investigate how the sensitivity can be improved by cross 
correlating two interferometers. The most general form for the cross correlation of 
two detectors is given by 

rT/2 rT/2 

C= dt dt' Si{t)S2{t')Q{t-t') 

J-T/2 J-T/2 

/oo roc _ _ 

df / df'6T{f - nsnminQif) . m 
-oo J —oo 

Here Q{t — t') is a filter function and the Sj's are the strain amplitudes that we read 
out from the i^^ detector: 

5,(t) = /if°(t) + s,,(t)+ni(t). (55) 
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The contributions are from resolvable astrophysical sources /if**™, the stochastic 
background Sj, and intrinsic detector noise rij. The function 5t that appears in the 
Fourier space version of the correlation function is the "finite time delta function" 

"'\te~--^^='^^. (56) 

-T/2 TIJ 



It obeys 

After removing resolvable astrophysical sources, the ensemble average of C is given 



5t(0) = T and 5{x) = lim 5r(a;) . (57) 

T^oo 



by 

{C) = {si, S2) + (si, ^2) + {rii, S2) + {ni,n2) = (si, S2) • (58) 

Notice that terms involving the noise vanish as the noise is uncorrelated with the signal, 
and the noise in different detectors is also uncorrelated. Thus, 

df / dfMf-ncswunQif')- (59) 

-00 J —00 

Our next task is to evaluate the quantity {s\{f)s2{f')) that appears in equation (|59|). 
Taking the Fourier transform of the expression in (p^ and performing the ensemble 
average we find 

{m)W)) = ls{f-f)S,{f)^{f)'^sm'(3, (60) 
where 7(/) is the overlap reduction function 

7(/) = J ^ {Fnn, /)F+(0, /) + rnn, m^Q, /)) e-2-/^ (---^)/^ . (61) 

The normalization is chosen so that 7(0) = 1. The function 7(/) is Hermitian in the 
sense that 712 (/) = 72i(/)- The inverse Fourier transform of 7(/) is a real function since 
7*(/) = 7(— /). The term "overlap reduction function" refers to the fact that 7(/) takes 
into account the misalignment and separation of the interferometers. In addition, 7(/) 
contains the transfer functions for each detector, as can be seen by setting Fi = F2 in 
(|6l|) and comparing with (^7|): 

5 

2sh?/3 

In other words, 7(/) is an overlap reduction function and detector transfer function all 
rolled into one. Inserting equation ( |60D into equation (|^) yields the expectation value 
of C: 

T roc _ 

{C) = -sm'p / df S,{fMf)Q{f) 

o J —00 

2T r°° 

= ,in'(3 df SMniif )Qif)] ■ (63) 
5 JO 

In the final line we have switched from the mathematically convenient range of 
integration / G (— c>o, 00) to the physically relevant range / G [0, 00). 



7i=2(/) = TT-T15 W)- (62) 
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The noise in a measurement of C is given by = C — (C), and the signal to noise 
ratio (squared) is given by 

Our task is to find the filter function Q{f) that maximizes this signal to noise ratio. A 
lengthy but straightforward calculation yields 

r2\ _ / I/O/' /M2 ; 



{N') = jJ_JQ{f)\'M{f)df. (65) 
where 

M(/) = 5„i(/) s„.2if) + s„,{f)s,{f)n2if) + ^„2(/)5,(/)7^l(/) 

Here Sn,i and TZi denote the spectral noise and the transfer function for the z*^ 
interferometer. The square of the signal to noise ratio can be written as 

SNR2 = — sin^/ji^^iSi^, (67) 
25 ^ {Q,Q} 



where {A, B} denotes the inner product | 

/CO 
dfA*{f)B{f)M{f), (68) 
-oo 

and 

The signal to noise ratio is maximized by choosing the optimal filter "parallel" to P. 
Since the normalization of Q drops out, we set Q{f) = -P(/)|lll, [H]- Using this filter, 
the optimal signal to noise ratio for the cross correlated interferometers is given by 

or, equivalent ly, as 

SNR= = H^Trd/Mgi(Z). (71) 



507r4 Jo ' pM{f) 
In the limit that the noise power in each interferometer is very much larger than the 
signal power we find 

507r4 Jo rSnl{f)Sn2{f) 

which recovers the expressions quoted by Flanagan and Allen |T5| when we set 
j3 = tt/2. The factor of sin"^/? was missed in these papers, as the implicit assumption 
that (3 = 7t/2 crept into supposedly general expressions. The main new ingredient in our 
expression are the transfer functions that reside in the overlap reduction function 
7(/). The transfer functions prevent us from performing the integral over the two-sphere 
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in equation (^Tj) in closed form except at low frequency. However, it is a simple matter 
to numerically evaluate 7(/) for a given detector configuration. 

If we are trying to detect a weak stochastic background with a noisy detector, then 
equation ([72| ) is the appropriate expression to use for the SNR. However, there may be 
some range of frequencies where the signal power dominates the noise power in each 
detector. The contribution to the SNR from a clean frequency window of this sort is 
given by 

SNR^(/, A/) = T ^ if 7 . . (73) 

25 i/-A//2 (h.(/0P + Rl(/0K2(/')/(2/5sm^,9)2 



The integrand is approximately equal to 1/2 for all frequencies, so that 



SNR(/,A/)^-sin2/5v'TA/. (74) 

This expression provides a useful lower bound for the SNR that can be achieved for 
detectors with a clean frequency window of width A/ centered at some frequency /. 
For example, the spectral density due to white dwarf binaries in our galaxy is expected 
to exceed LISA's noise spectral density for frequencies in the range 10"'^ 3 x 10^^ 
Hz. To detect this background at 90% confidence requires a signal to noise ratio of 
SNR=1.65, which can be achieved with less than three hours of integration time T. 



J^..^. Cross correlating two LISA interferometers 

One possible modification to the current LISA proposal ^ would be to fly six spacecraft 
instead of three, and thereby form two independent interferometers. The cost of doing 
this is considerably less than twice the cost of the current proposal, and it has the 
advantage of providing additional redundancy to the mission. How would we best 
use a pair of LISA interferometers? If our main concern is getting better positional 
information on bright astrophysical sources, then we would fly the two interferometers 
far apart, eg. with one leading and the other trailing the Earth. However, if we want 
to maximize the cross correlation then we need the interferometers to be coincident 
and coaligned. However, a configuration of this type is likely to share correlated noise 
sources, which would defeat the purpose of cross correlating the interferometers. A 
better choice is to use a configuration that is coaligned but not coincident. This can be 
done by placing six spacecraft at the corners of a regular hexagon as shown in Figure 
5. Notice that interferometers have parallel arms and that the corner spacecraft are 
separated by the diameter of the circle. 

Using the same coordinate system that we used earlier for a single interferometer, 
the unit vectors Uj and Vj along each interferometer arm are: 

1 . v^. 

ui = -U2 = 2^ + -^y 

vi = -V2 = + — y, (75) 
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Figure 5. The hexagonal cross correlation pattern. 



where the index i = 1,2 labels the two interferometers. The displacement of the corner 
spacecraft is given by 

xi-X2 = -2i?y (76) 

where R is the radius of the orbit. The angle between the interferometer arms is 
(3 = 7r/3, and the length of each arm is L = y/3R. The light crossing time between 
the two interferometers, 2R/c, is almost equal to the light crossing time along each 
interferometer arm, v^-R/c. Thus, the loss of sensitivity due to multiple wavelengths 
fitting between the interferometers occurs for frequencies near the transfer frequency 
f^: (the transfer frequency corresponds to wavelengths that fit inside the interferometer 
arms). For LISA the transfer frequency is /* = 9.54 x 10^^ Hz. 



The various ingredients we need to calculate 7(/) follow from those given in 
equations ( pO]) through (0). Putting everything together in (|6ll) and working in the 
low frequency limit we find 



143700480 V/, 

For frequencies above the overlap reduction function decays as /~^. A numerically 
generated plot of 7(/) is displayed in Figure 6. Scaling the signal to noise ratio (|72p in 
units appropriate to a pair of LISA interferometers we have 

SNR^ . 34,4 I ^] r (^] h(/)P ('-^J (^J (^)' , (78) 




Using equation ( pi]) for S'„(/) and assuming a scale invariant stochastic background 
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Figure 6. The overlap reduction function for the hexagonal cross correlation 
patern. The solid line was generated numerically while the dashed line is our analytic 
approximation from equation j?^). 

yields a signal to noise of 

SNR = 1.44./^ { . (79) 
V year V lO-^^ j ^ ' 

This represents a 500 fold improvement on the sensitivity of a single LISA detector. 
Indeed, if it were not for the astrophysical foregrounds, a pair of LISA detectors would 
be well poised to detect a scale invariant gravitational wave background from inflation. 

5. Missions to detect the CGB 

We will work on the assumption that astrophysical foregrounds swamp the CGB for 
frequencies above a few /zHz, and design our missions accordingly. What we have in mind 
is a post-LISA mission based on the (by then) tried and tested LISA technology, but with 
some allowance for improvements in basic technologies such as the accelerometers. Since 
cross-correlating two interferometers at ultra-low frequencies does not buy us a major 
improvement in sensitivity, we need to start with a design that has excellent sensitivity 
at low frequencies. Basically this means building bigger interferometers with better 
accelerometers. To be concrete, we show the sensitivity curves for three generations of 
LISA missions in Figure 7. LISA I corresponds to the current LISA design with the 
spacecraft cart-wheeling about the Sun at 1 AU, separated by L = 5 x 10^ m. LISA II 
refers to a possible follow on mission with the three spacecraft evenly spaced around an 
orbit at 1 AU, so that the spacecraft are separated by L = -\/3 AU. The LISA II mission 
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would use similar optics||to LISA I (same laser power and telescope size), but allows for 
an order of magnitude improvement in accelerometer performance. LISA III is similar 
to LISA II, except that the constellation would orbit at 35 AU (between Neptune and 
Pluto) and the accelerometers would be improved by a further two orders of magnitude. 
The acceleration noise for LISA III would benefit from the three orders of magnitude 
reduction in solar radiation relative to LISAs I and II, but this would come at the cost 
of having to power the spacecraft using nuclear generators (RTGs). 




For all three missions we see that the /iHz range lies well below the interferometer's 
transfer frequency 

f, = lMxlQ-'[^^^ Hz. (80) 

Indeed, it would take a mission orbiting at 180 AU ( somewhere inside the Kuiper belt) 
to achieve a transfer frequency of = 10~^ Hz. In practical terms this allows us to 
work in the low frequency limit / ^ /* where it is easy to derive analytic expressions 
for the overlap reduction function 7(/). Moreover, we need only consider acceleration 
noise when working below 1/iHz. All our calculations are based on the same hexagonal 

The main difference is that we have to "lead our target" by a much larger amount for the LISA 
II mission. In other words, the angle, A0, between the received and transmitted laser beams is 
much larger for LISA II than for LISA I. A simple calculation yields /S.9 — 2w/c, where v is the 
velocity around the circle shown in Figure 5. For LISA II v/c — [GMqI c? j K)^^"^ , while for LISA I 
vjc = 2e{GMQ/c^ /R)'^/'^. Here e = 0.01 is the eccentricity of the LISA I orbits and i? = 1 AU. This 
equates to lead of = 4 x lO"*^ radians for LISA I and = 2 x lO"'' radians for LISA II. 
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cross correlation shown in Figure 5 that we used for the LISA I cross correlation. Using 
equation (|50| ) for the acceleration noise and ignoring the position noise yields 

In the limit that the noise power dominates the spectral density of the CGB we arrive 
at an estimate for the signal to noise ratio for frequencies below l/iHz 

SNR^ = 29.4 I — I I ^ r X 

The SNR is scaled against the LISA II specifications. Since ///* ^ 1 it is a good 
approximation to set 7(/) = 1 inside the integral. We model the CGB spectrum near 
1/iHz by a simple power law: 

^^gw(/) = ^^^Hz , (83) 

scaled relative to a reference value at l/xHz. With these approximations we arrive at our 
final expression for the SNR for a pair of cross correlated interferometers in the ultra 
low frequency regime: 

3.13 / T L \ V3 X 10-l6ms-2^' /n ;,2\ / , \ "+3/2 




SNR 




We can compare this to the SNR of a single interferometer in the ultra low frequency 
regime by combining equations (PBT) and (F^) to find 




SNR.(/) ^ 0,683 - M I I ^ 1 . (85) 

We see that the individual and cross correlated sensitivities are related: 

SNR = , (TA/)^/2 SNR,(/) , for / ^ A/ ^ 1 /iHz . (86) 

'I + 2a/3 



This supports our earlier assertion that cross correlating two interferometers improves 
on the sensitivity of a single interferometer by (TA/)^/^ where A/ is approximately 
equal to the central observing frequency /. 

Equation ( p4D tells us that a GABI detector built from two LISA II interferometers 
could detect the CGB at greater than 90% confidence with one year of data taking if 
f2gw(l/xHz)/zQ = lO"^''. The same equation also tells us what has to be done to achieve 
greater sensitivity. It is clear that increasing the duration of the mission is not the best 
answer as the SNR only improves as the square root of the observation time. In contrast, 
increasing the size of the interferometer or reducing the acceleration noise produces a 
quadratic increase in sensitivity. For example, a GABI detector built from a pair of 
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LISA III interferometers could detect the CGB at 90% confidence for signals as small 
as ng„(l^Hz)/iQ = 4 X 10~^^. In-fact, for signal strengths above this value the LISA 
III detectors are signal dominated and we need to use equation ( ffil) instead of ( PH) to 
calculate the SNR. 

We can also use equation (Q) to determine how well a GABI mission can measure 
the CGB power spectrum. Suppose that we break up the frequency spectrum into bins 
of width 5f. Nyquist's theorem tells us that 5f > 1/T, eg. for an observation time of 
one year the frequency resolution is 3.17 x 10^*^ Hz. Writing 6f = n/T where n > 1, 
and taking 6f ^ /, we find from (|8^) that in the frequency window (/ — 5//2, f + 6f /2) 
the SNR is 

We see that the SNR in each frequency bin scales as Z^"*"". For a scale invariant spectrum 
{a = 0) this translates into poor performance at frequencies below 1 /xHz and limits 
the range over which we can measure the spectrum. For example, if the CGB has a 
scale invariant spectrum and an amplitude of ^Igwh^ > 2 x 10~^^, a pair of LISA II 
interferometers could measure the spectrum over the range 10~^ — > 10~^ Hz. For a 
pair of LISA HI detectors the main limitation at low frequencies comes from Nyquist's 
theorem. Since mission lifetimes are limited to tens of years, it will not be possible 
to measure the CGB spectrum much below 10^*^ Hz using space based interferometers. 
Indeed, it may be difficult to push much below 10^^ Hz unless ways can be found to 
build detectors that are stable for many months. A pair of LISA HI interferometers 
could measure the spectrum between 10~^ 10"^ Hz if i7g„(10^^ Hz)/;,q > 10~^^. To 
measure the CGB spectrum below 10"*^ Hz requires a return, full circle, to the world of 



CMB physics. Detailed polarization measurements of the CMB can be used to infer 
the CGB power spectrum for frequencies in the range 10"^*^ —>■ 10~^^ Hz. 



6. Discussion 



The LISA follow-on missions we have described will be able to detect or place stringent 
bounds on the CGB amplitude and spectrum between 10~^ and 10~^ Hz. But how can 
we be sure that it is the CGB we have detected and not some unresolved astrophysical 
foreground? The answer can be found in the statistical character of the competing 
signals. Most early universe theories predict that the CGB is truly stochastic. In 
contrast, the astrophysical signal is only approximately stochastic, in a sense that can be 
made precise by appealing to the central limit theorem. The line in Figure 1 indicating 
the predicted amplitude of various astrophysical foregrounds corresponds to what is 
known as the confusion limit. It marks the amplitude at which we can expect to find, 
on average, one source per frequency bin[p|"|. It is not until we reach signal strengths 

ft The confusion limit goes down when the frequency resolution goes up. Most plots of the confusion 
limit assume a one year observation period so the bins are 3.17 x 10~* Hz in width. 
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considerably smaller than the confusion limit that the astrophysical signal starts to look 
stochastic. And therein lies the answer to our question. A "stochastic" signal due to 
astrophysical sources will always have bright outliers that sit just above the amplitude 
of the stochastic signal. The number and distribution of the outliers can be predicted 
on statistical grounds. If we do not see the outliers, then we can safely conclude that 
we have detected the CGB. 
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